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The importance of nuclear reactions in low-density astrophysical plasmas with ion tem- 
peratures T > 10 10 K has been recognized for more than thirty years. However, the lack 
of comprehensive data banks of relevant nuclear reactions and the limited computational 
power have not previously allowed detailed theoretical studies. Recent developments in 
these areas make it timely to conduct comprehensive studies on the nuclear properties 
of very hot plasmas formed around compact relativistic objects such as black holes and 
neutron stars. Such studies are of great interest in the context of scientific programs of 
future low-energy cosmic 7-ray spectrometry. In this work, using the publicly available 
code TALYS, we have built a large nuclear network relevant for temperatures exceeding 
10 10 K. We have studied the evolution of the chemical composition and accompanying 
prompt gamma-ray emission of such high temperature plasmas. We present the results 
on the abundances of light elements D, T, 3 He, 4 He, 6 Li, 7 Li 9 Be, 10 B, 11 B, and briefly 
discuss their implications on the astrophysical abundances of these elements. 
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1. Introduction 



Low-density optically thin, hot plasmas with ion temperature exceeding 1 MeV 
can form near compact relativistic objects, such as in accretion flows close to black 
holeJ^and in strong shock waves related, for example, to supernova explosions PI The 
most straightforward approach to studying the dynamics of formation and evolution 
of such plasmas is the detection of prompt 7-ray lines which are unambiguous 
signatures of specific nuclei. 

Collisions in thermal plasmas with temperatures exceeding 1 MeV are character- 
ized not only by excitation reactions with emission of gamma-ray lines, but also by 
the destruction of nuclei. Generally, in such plasmas, the destruction of nuclei pro- 
ceeds on timescales shorter than excitations ones. Thus, only prompt de-excitation 
lines related to specific nuclei are expected to contribute to the plasma radiationP 
The presence of a given prompt 7-ray line is related to the existence of a given nu- 
cleus, which is itself governed by the destruction processes. This suggests that there 
should exist a strong relation between the gamma-ray line spectrum and nuclear 
destruction processes. 

The break-up processes are sensitive to the plasma temperature and density. 
Therefore, the detection of 7-ray lines and their spectral profile evolution may reveal 
unique information about the chemical composition and physical parameters of the 
plasma, such as temperature and density at any given time. 

The first attempts to assess the importance of nuclear processes in hot, thin 
astrophysical plasmas were made in Refs. |3j |4j Different astrophysical implications 
of nuclear reactions in hot two-temperature plasmas have been discussed also in 
Ref. [5] - [131 These works provide estimates for the isotope ratios of different light 
elements such as Deuterium (D), 6 Li, 7 Li 9 Be, 11 B, etc, and also estimates for 
the 7-ray line luminosity due to nuclear reactions. A mechanism to produce light 
elements, alternative to Big-Bang Nucleosynthesis, is also suggested in Refs. [3l [4] 
In particular, they discuss a mechanism for deuterium production in hot plasmas. 

In these early works, however, only a limited number of reactions could be 
taken into account, mainly due to the lack of data on nuclear cross sections and the 
increased computational complexity of calculations as more nuclei are considered. 

Most of these difficulties have nowadays been overcome. The introduction of 
on-line nuclear database d^^ provides easy access to experimental and theoretical 
data on nuclear cross sections. There also exist publicly available computer codes 
that allow the calculation of nuclear cross sections and other quantities of interest. 
These make use of several theoretical models. 

In this work we build a large nuclear reaction network and use it to calculate the 
temporal evolution of the composition, the 7-ray line emissivity, and the correspond- 
ing 7-ray spectrum, for a low density plasma with ion temperature T > 1 MeV. 
Two cases are considered: when the temperature and density evolve with time, 
and when they remain constant. We present results for the abundance of light 
elements and the baryonic 7-ray spectrum due to nuclear reactions and neutron- 
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proton Bremsstrahlvmg. Finally, we discuss some mechanisms that could enhance 
the production of light elements. This could be of special importance, because it 
could modify the cosmic abundance of light elements predicted by the Standard 
Big-Bang Nucleosynthesis model. The production of neutrons is also of consider- 
able interest, since these particles can escape from the source and interact with 
nearby targets. 

2. Cross sections and reaction rates 

The calculation of the evolution of the nuclear abundances and the 7-ray line emis- 
sion in low density, high temperature plasmas, requires the use of a comprehensive 
nuclear network. This network should take into account all relevant nuclear reactions 
that may occur in a hot plasma with T > 1 MeV. 

There is an extensive literature about how to build and solve nuclear networks, 
e.g. see Refs. [THI - 0£H Nowadays, all this knowledge is compressed and very easy 
to access through on-line nuclear databases for astrophysical applications, such as 
JINA Reaclib^ NON-SMOKER^ NucAstroDataP® etc. There also exist many 

compute r pr ograms which store, manage and solve the networks, e.g. libnuceq and 
1071 

libnucnet¥ LLi However, all these nuclear networks are designed to be applied in sce- 
narios such as big-bang nucleosynthesis, core collapse supernova, stellar interiors, 
etc. Generally, they are valid for temperatures T < 10 10 K. Therefore, to describe 
nuclear reactions for astrophysical plasmas with T > 10 10 K (or T > 1 MeV), we 
find it necessary to extend the aforementioned works. 

Without loosing generality, our main purpose is to describe the evolution of the 
abundance of nuclear species and the 7-ray emission of a very hot, low density and 
optically thin astrophysical plasma. To further simplify the problem, we assume the 
following conditions hold: 

• all nuclei are instantaneously thermalized, and their velocity distributions 
are described by a Maxwellian distribution function at temperature T, 

• all target nuclei interact in their ground state, 

• the plasma has low density (n ~ 10 18 cm -3 ) and high temperature (T > 
1 MeV), and its initial composition is similar to the standard cosmic/solar 
chemical abundances. 

These assumptions dramatically simplify the reaction network. The first assump- 
tion reduces the calculation of the reaction rates to a simply Maxwellian average 
< av >. Thus, < av > will be only a function of the temperature T as shown in 
Ref.HU 

< " ^ (T) " fjj^f f «» (-isO dE - (1) 

where E is the projectile CM-frame energy, cr^ is the cross section of the binary 
interaction i + j — > k + jiij is the reduced mass of the interacting system, and 
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is the Boltzmann constant. 

The second assumption simplifies both the calculation of the cross sections and 
the reaction rates. In a high temperature plasma, different isomeric states of a given 
nucleus may be in thermal equilibrium. There fore , to obtain the reaction rate, one 
should sum the contributions of all isomers 112II2QI-|30| ^ s a £ rg ^ a pp r0 ach, we have 
assumed here that all isomeric states decay before the nuclei interact. 

The third assumption restricts the number of targets/channels and interactions 
to be taken into account. 

Under these assumptions, and since we are interested in low density plasmas, 
only the binary interactions contribute to the chain of nuclear reactions. For the 
composition of plasmas similar to standard astrophysical environments, only the 
reactions with participation of protons (p), neutrons (n), and alpha particles (a) 
are important, c.f. |Appcndix A.l| 

Although the nuclear reaction rates for p, n and a projectiles can be simply 
calculated from the local conditions of the plasma - such as ion temperature and 
density - this is not the case for the interactions of photons with nuclei. The 7- 
interaction reaction rates not only depend on local parameters, but also on the 
number density of photons, which itself depends on the optical depth of the plasma. 
Furthermore, the optical depth is also a function of "non-local" variables such as 
the electron number density and size of the plasma. 

The 7-rays are produced and transported inside the plasma. Hot two- 
temperature plasmas radiate X-rays and 7-rays through interactions of both com- 
ponents of the plasma. The dominant radiative cooling channels for ions are the 
excitation of nuclei with production of prompt 7-ray lines, the 7-radiation released 
at p — n capture, as well as p — n Bremsstrahlung. The electronic component radiates 
more efficiently, predominantly through Bremsstrahlung, and in the case of large 
optical depth, is cooled due to comptonization. 

Because the photon number density depends on the optical depth of the plasma, 
the effect of the photodisintegration of nuclei cannot be estimated without consid- 
ering the plasma as a whole. The importance of this effect must be checked for each 
specific problem. In this paper we have neglected photodisintegration. Generally, 
this is a good approximation for optically thin, low density plasmas. We illustrate 
in |Appendix A. 2 the effect of the photodisintegration on 16 and Deuterium nu- 



clei, for the case of a spherically symmetric and optically thin plasma. We show that 
even if the electron temperature is comparable with the ion temperature, the effect 
of photodisintegration on the evolution of the plasma can still be safely neglected. 
However, this might not be the case for an optically thick plasma, or for the case 
where nuclear statistical equilibrium is reached. 

The integral in Eq.(fT|) is evaluated over the energy range (0, 00). In practice, 
however, one needs to integrate Eq.(fT|) only in a "narrow" energy interval from which 
comes the most important contribution to the reaction rate. To have an accurate 
estimation of the reaction rate we must know precisely the cross sections in that 
narrow energy range, since the contribution of the leftover interval is negligible, by 
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definition. As discussed in | Appendix A. 3 in the temperature range 0.1-20 MeV the 



most relevant interval of projectile energy is that between 0.01 MeV and 150 MeV. 

The Maxwellian average < a v > is the main ingredient to solve the kinetic 
equations It is convenient though, to work with mass fraction abundances instead 
of number densities. If we define the mass fraction abundance of the element i as 
X 1 — Ai n % I p, where A4 is the mass number, and p = A\ n l the nucleon density^ 
then for our case the kinetic equation reads: 

n,p,a N / a \ 



The index i runs over the three types of projectiles (protons, neutrons and alpha 
particles), whereas j runs over all TV elements of the network. 

The remaining step for calculating < a v >y in Eq.([T]), is to provide the cross sec- 
tions in the rele vant projectile energy interval. We chose to use the public available 
code TALYSpS] to calculate the cross sections tr% for all relevant nuclear interac- 
tions, together with the 7-ray lines associated with these interactions. TALYS is an 
attractive tool to be applied to the study of nuclear processes in high-temperature 
plasmas because: 

i. the energy of the projectile can be in the range 1 keV-200 MeV, 

ii. it incorporates modern nuclear models for the optical model, level densities, 

direct reactions, compound reactions, pre-equilibrium reactions, 

iii. it can treat 7, n, p, D, T, 3 He, and a-particles as both projectiles and ejectiles, 
vi. target mass numbers can be between A = 12 and A = 339. The code can also 

calculate cross section for targets with mass number in the range A = 6 — 12, 
but the results may not be as accurate. 

We used TALYS to calculate the cross sections for reactions involving elements 
from 6 Li to 70 Zn. For targets with A < 6 (i.e. n, p, D, T, 3 He, 4 He), we collected 
the cross sections from theoretical evaluations; when no evaluations were available 
we used experimental data. Both evaluations an d ex perimental data points were 
taken from the publicly available databases ENDpISl and EXFORp^ 

The reaction network should in principle satisfy a "closure condition" : the reac- 
tion products of all nuclei in the network must be also part of the same network. 
In practice we select the most relevant reactions/channels, which implies that the 
closure condition holds only approximately. We can fix the accuracy up to which the 
closure condition is satisfied by introducing a threshold value for the Maxwellian 
average. Every channel with a Maxwellian average smaller than the threshold is ig- 
nored. Changing the value of the threshold changes the number of reaction channels 
in consideration, therefore this parameter plays a crucial role in the accuracy of the 
solution, the calculation time, the dimension of the network, etc. 



a In practice, nucleon density p can be converted to mass density pmass , simply by multiplying it 
with nucleon 
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To calculate the cross sections and reaction rates for nuclei from 6 Li to 70 Zn 
we proceeded as follows. Using TALYS we calculated the cross sections for the 
interaction of solar abundance nuclei (the stable nuclei) from 6 Li to 70 Zn, with 
n, p and a as projectiles. Then, the corresponding reaction rates were calculated 
as in Eq.(TT|). At T = 1 MeV the largest reaction rates are of the order of Rfj ~ 
10~ 16 — 10~ 15 cm 3 s~ 1 . We chose a threshold value R th = 10 -19 cm 3 s _1 , ensuring 
that the contribution of the dismissed reactions would be no more than 0.1%. 

From the selected channels new targets may appear. They are potentially created 
in the plasma due to the reactions of the previous generation of nuclei, and enter 
in the calculations as a second target generation. We then repeated the previous 
procedure for the new targets, i.e. calculating the cross sections, the reaction rates, 
and the selection of the relevant channels according to the chosen threshold value. 
The new channels may still introduce a third generation of new nuclei. We applied 
the algorithm until there is no generation of new nuclear species. The network is 
then closed up to the accuracy given by R th . Notice that the value of R th must be 
adjusted according to the scope of the problem, and there is no general advice for 
its choice. 

We extended our calculations by fitting the Maxwellian averages < av >^ as a 
function of temperature using the parameterization with seven parameters suggested 
in Ref . H3 Ref. M ■ 



<av >* = exp(ai + a 2 T" 1 + a 3 + ai T^ + a 5 T + a 6 T^ + a 7 log T). (3) 



The fit is valid for T £ [0.1, 20] MeV. We estimated the goodness of the fit for the 
seven parameters using the x 2 method. In the mentioned temperature range the 
majority of the rates have a residual fluctuation less than 5%. 

As we have already mentioned, for light elements with A < 6 we obtained the 
relevant cross sections from evaluations or, in the case where no evaluations were 
found, from experimental data points. It must be mentioned here that using eval- 
uations or experimental data introduces additional uncertainties. Some evaluations 
do not exceed projectile energies of 20 MeV. In these cases, we extrapolated the 
cross section in the high energy region as a ~ 1/E, see e.g. Ref. HH Listed below 
are all the reactions involving nuclei with A < 6 we have included in this work: 
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HI) X H (n, 7) D 

Dl) D (n, 2n) X H 
D2) D (n, 7) T 

Tl) T (n, 2n) D 
T2) T (n, 3n) 4 H 

hi) 3 He (n, p) T 
h2) 3 He (n, D) D 

al) 4 Hc (n, D) T 
a2) 4 He (p, D) 3 Hc 



D3) D (p, np) : H 
D4) D (p, 7) 3 He 

T3) T (p, n) 3 He 
T4) T (p, D) D 

h3) 3 He (n, 7) 4 Hc 
h4) 3 He (p, 2p) D 

a4) 4 He (a, n) 7 Bc 
a3) 4 He (a, p) 7 Li 



D5) D (a, np) 4 He 



T5) T (p, 7) 4 He 
T6) T (a, n) 6 Li 

h5) 3 He (a, 7) 7 Be 



a5) 4 He (a, 2p) 6 He 
a6) 4 He (a, np) 6 Li 
a7) 4 He (a, D) 6 Li 



The data on the corresponding cross sections were taken from the following sources: 

HI) Reaction cross-section from ENDF/B-VI.8 evaluations H^l valid for a pro- 
jectile energy up to 150 MeV. 
D1-D5) Data from ENDF/B-VII.f|IS valid for projectile energy up to 150 MeV. 
T1,T3,T4) Data from ENDF/B-VII.l evaluations^ valid for projectile energy up to 
20 MeV. 

T2) Data from JEFF-3.1 evaluations^!, valid for projectile energy up to 20 
MeV. 

T5) For this reaction we have not found any theoretical evaluation, therefore 
we have used the experimental dataP^From 10-100 keV in CM-system, we 
have taken the data from Canon et. al. (2002)P1 For 8.34 and 13.6 MeV 
we got it from Calarco et. al. (1983) 

T6) Again no evaluations, we have extracted the data from R. Spiger and T. 
Tombrello (1967)SH These data are valid up to 18 MeV. 
hl-h4) Evaluations from ENDF/B-VII.O, which are valid up to 20 MeV. 

h5) No evaluations. Experimental data in the interval 0.02-3.13 MeV collected 
by different authors (experiments from 2008-2009). 

al) Data from R. Shamu and J. Jenkin 1964 around 22Mev32 and one point 
at 50MeV differential cross sections from Sagle et. al. (1991)031 

a.2) Experimental data for 31 MeV from Bunch et. al. 1964,351 for 32 MeV, 40 
MeV, 50 MeV and 52.5 MeV from Sagle et. al. (1991) S3 for 85 MeV from 
Votta et. al. (1974)33 and for 200 MeV from Alons et. al. (1986)33 
a3-a7) Data taken from King et. al. (1977) for energies between threshold and 50 
M e yu±oi p or ener gj eg higher than 60 MeV, measurements from Glagola et. 
al. (1982)321 and fits from Mercer et. al. (2001) were usedP^In the work of 
King et. al. (1977), there are also data for the excited-state cross section of 
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7 Li(0.478) production 7-ray line. 

Finally, it is important to note that nuclear decay becomes a non-negligible effect 
over sufficiently long timescales. In particular, neutron decay plays a major rol e in 
our problem due to its high abundance and relatively short lifetime, t„ — 885.7 sP^l 
We have taken into account neutron decay in our calculations. 

3. Network solutions 

After calculating all the < av we are ready to solve the system of differential 
equations Eq.Q for each species in the network. Our network contains a relatively 
large number of nuclear species (around 270). Therefore the system of differential 
equations is too complicated to be solved analytically, and must be solved numer- 
ically. We used the Matlab®'s "odel5s" solver, which is an implementation of the 
Klopfenstein-Shampine family of Numerical Differentiation Formulas (NDF) of or- 
ders 1 to 5.^ 

The numerical solver does not guarantee a physically meaningful solution, we 
then included some constraints to be checked in every solution step. The constraints 
originate from conserved quantities such as the number of nucleons (or baryon 
number conservation), electric charge, etc. One of the most important constraints 
is baryon number conservation, which can be written as 

P = => J2*i = ° ^Aihi = Q, (4) 

where p = ^ AiJii is the number density of nucleons. This condition ensures that 
the number of nucleons remains constant even if, for instance, neutrons decay. 

3.1. Evolution of the plasma chemical abundance and the j-ray 
emission 

We present the solution of the nuclear network for some simple cases. As a first 
example, we considered a unit volume of plasma at two fixed-temperatures, 1 MeV 
and 10 MeV. The nucleon number density for this plasma is p — 10 18 cm~ 3 . We 
studied three special cases, each of them differing only in the initial plasma compo- 
sition. The initial mass abundance compositions for each case are: 

1) X p = 0.7 and Xi2 C = 0.3 

2) X p = 0.7 and A 2 8 Si = 0.3 

3) X p = 0.7 and X56 Fc = 0.3 

This is equivalent to having plasma with n p = 7xl0 17 cm -3 and ni2 C /n p = 3.6 10~ 2 , 
n2s Si /n p — 1.53 10 -2 , and n56p c /n p = 7.7 10~ 3 , respectively. We then evolved the 
plasma for a period of 100 s. The temporal evolution of the plasma composition and 
the 7-ray line emission is shown in Figs. [T] to [3l Further discussion on these results 
is presented in Section 01 
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Fig. 1. Temporal evolution of mass-abundances and 7-ray line emission for a plasma with initial 
composition of 70% protons and 30% 12 C, for temperatures 1 MeV (left) and 10 MeV (right). In 
each case, the bottom panels show in particular the evolution of the abundance of light elements 
up to 11 B, multiplied by the coefficient in brackets for comparison. 
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Fig. 2. Temporal evolution of mass-abundances and 7-ray line emission for a plasma with initial 
composition of 70% protons and 30% 28 Si, for temperatures 1 Mev (left) and 10 MeV (right). In 
each case, the bottom panels show in particular the evolution of the abundance of light elements 
up to 11 B, multiplied by the coefficient in brackets for comparison. 
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Fig. 3. Temporal evolution of the mass-abundances and 7-ray line emission for a plasma with 
initial composition 70% protons and 30% 56 Fe, for temperatures 1 Mev (left) and 10 MeV (right). 
In each case, the bottom panels show in particular the evolution of the abundance of light elements 
up to 11 B, multiplied by the coefficient in brackets for comparison. 
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3.2. The plasma's baryonic component -f-ray spectra 

We calculated the 7-ray spectrum result only of baryonic interactions. For high 
temperature and low density plasmas, the dominant baryonic 7-ray emission pro- 
cesses are nuclear reactions and proton- neutron (p-n) Bremsstrahlung, see Ref.(j3]). 
Proton-neutron Bremsstrahlung and nuclear reactions such are the capture pro- 
cesses 1 H(n,7)D, D(n,7)T, D(jj, 7) 3 He, and T(p, 7) 4 He, give rise to the contin- 
uum component of the spectrum. The discrete component is purely due to the de- 
excitation of nuclei, that decay into different isomeric states leading to the emission 
of specific 7-ray lines. 

We define $ 7 as the gamma-ray emissivity per unit energy (or unit frequency) . 
If h-y is the number of 7-rays emitted per unit time and volume at energy E^ , then 
the emissivity can be written as follows: 



$ 7 (£ 7 ) = £ n 



dE~ 



Gamma-rays are directly created in capture reactions of the form 1 + 2 
such &s n + p — ;» D + 7. The emissivity for a capture reaction is given by 



(5) 
3 + 7, 



where 



<& 7 (.E 7 ) = n\ n% c 



m 3 



TTli + 7712 1/ 7T fl C 2 (fcfiT) 



2 H-,,T\3 E "i 



(6) 



I{E 7 ) = / dE 



E(E + Q)a(E) 



exp 



-.2„2 



2(mi + TO 2 )fc B T 



1 - 



E + Q 
E~. 



E 



Here ni, ni are the number densities of species 1 and 2, fi is the reduced mass of 
the interacting system (nuclei 1 and 2), £? 7 is the 7-ray energy in the LAB-frame, 
E is the collision energy in the center-of-mass frame, and cr(E) and Q are the cross 
section and the Q- value of the reaction, respectively. 

The recoil of nuclei 3 is neglected in Eq.©, i.e. terms oc E^/2m 3 are not in- 
cluded. The energy of the interacting particles is ~ feT, so most of the pho- 
tons will be produced with an energy of the order of E 1 ~ Q + ks T. Thus 
(Q + ks T)/ni3 < 10~ 2 and the approximation is justified. 

The discrete gamma-ray spectrum is, as we have already mentioned, the result 
of the de-excitation of nuclei. The excited nuclei have a Maxwellian velocity distri- 
bution, therefore the emitted 7-ray lines will be Doppler-broadened. The Doppler 
broadening of a 7-ray line at energy E® emitted by a nuclei of mass m in a plasma 
at temperature T is 



81n(2)fc B T 

— — ^7- 



(7) 
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In a high-temperature plasma is much larger than the natural width Tq for 
almost all the strongest lines. This simplifies the broadening to a simple Gaussian 
profile. If the energy of the photon emitted by the nucleus at rest is E®, and if n 7 is 
the number of photons with energy E® that are emitted per unit time and volume, 
then the emissivity can be written as 



d> 7 (£ 7 ) 



exp 



(E~ t -E°f 
2a% 



(8) 



where the Gaussian broadening <jq is 



lk B T 

° G = \ 2 E ~f 

V mr ' 

Finally, for the calculation of the emissivity due to n-p Bremsstrahlung we em- 
ployed the cross section together with the calculations from Ref.©. 

We now show the spectra on three specific moments during the plasma evolution, 
for the 12 C, 28 Si and 56 Fe solutions that we have already presented in Section [3l 
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Fig. 4. Gamma-ray spectra of a plasma with initial composition of 70% protons and 30% 12 C, 
for temperatures 1 MeV (left) and 10 MeV (right). 



3.3. Other interesting examples 

We now show, for comparison, two other interesting examples. The first corresponds 
to the case when the plasma temperature is assumed to be constant over the entire 
solution interval; in the second case the plasma cools on timescales comparable to 
nuclear destruction ones. 

To a first-order approximation, the result of the evolution of a low-density, high- 
temperature plasma over timescales larger than one second is a system formed by 
protons and neutrons. Since free neutrons decay, the plasma eventually consists 
purely of protons. 

The other case considered here is that where the plasma cools "instantly": if 
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Fig. 5. Gamma-ray spectra for the plasma with initial composition of 70% protons and 30% 28 Si, 
for temperatures 1 MeV (left) and 10 MeV (right). 

ta is the typical cooling time and r e the evolution timescale, then ta <r e . Under 
this conditions many light and intermediate-mass elements can form. There is no 
physical impetus for this example; rather we performed this calculation to illustrate 
that the freezing out of the reactions with a change in temperature does indeed 
occur. However, a physically plausible scenario for such a change in temperature 
could be a plasma interacting with a fast shock. The same phenomenon could occur 
if the density and/or pressure are instantaneously lowered, such as in the case of 
free expansion. 

We adopted a temporal dependence for temperature of the form T = 0.1 + 
20 exp(— i/0.001) MeV to illustrate the freezing of reactions. Due to the added 
complexity and computational expenses, our network computes reaction rates down 
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to T = 0.1 MeV, and therefore we do not observe the total freezing. 
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Fig. 7. The temporal evolution of the mass fraction abundances and the 7-ray line, emission for 
a plasma with initial composition of 70% protons and 30% 56 Fe, for constant temperature 10 MeV 
(left), and for time dependent temperature T = 0.1 + 20 exp(-i/0.001) MeV (right). 



4. Discussion 

As discussed above, we validated our network by calculating the mass abundance 
evolution of all the nuclei for some special cases. Figures Q] to [3] show the temporal 
evolution of the mass fraction abundances and the number of photons of a given 
energy emitted per unit time and volume for the cases considered. Also shown is 
the number density ratio with respect to protons (ni/n p ) for the light elements n, 
D, T, 3 He, 4 He and 6 Li, 7 Li, 9 Be, 10 Be, and 10 B, U B. 

By comparing the plots on the left and on the right, we can draw some general 
inferences. It is clear that the value of the temperature has a strong impact on the 
characteristic timescales of the reactions. This is because the rates of almost all 
reactions increase with temperature, and correspondingly the nuclear abundances 
and the emission of 7-rays evolves much faster. Note from Figures Q] to [3] that some 
intermediate elements (the result of break-up processes) , become more abundant in 
certain periods of the evolution. At higher temperatures more interaction channels 
are opened, which leads to more 7-ray lines due to the de-excitation of nuclei. In the 
case when the plasma cooling timescale is such that all reactions freeze at a time 
when intermediate elements are relativity more abundant, this can result in higher- 
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than-normal abundances of such elements at the end of the network's evolution, as 
shown in Fig. [7J 

The tendency of T > 1 McV plasmas is to destroy, very quickly, all nuclei 
heavier than 4 He, and to a first approximation, we may conclude that on time 
scales r ~ 1 s a proton-neutron plasma is formed. We also observe for times t > 1 s 
that some small traces of deuterium, tritium, 3 He and 4 He are present. This is a 
result of the break-up processes and weak fusion reactions which take place on longer 
timescales. A result of the proton-neutron plasma formation is that all 7-ray line 
emission disappears. For timescales r > 1 s, neutron capture by protons through 
the reaction 1 H(n, 7)D, can produce deuterium and the emission of 2.22 MeV 7-ray 
photons. Due to the kinematics the 2.22 MeV emission line should be broadened 
and shifted - by approximately ksT - as is shown in Figs. 0] to Fig. [6l 

Figures QJ3] also show the temporal evolution of the number density ratio, nii/n p , 
of stable light elements with mass numbers up to A = 11. The abundance of light 
elements such as Li, Be and B is practically negligible for a constant temperature 
and density plasma. Hydrogen and helium isotopes, however, are much more abun- 
dant due to proton and neutron equilibrium, provided that the system remains in 
equilibrium. 

High temperature, low density plasmas have a characteristic 7-ray spectrum of 
baryonic origin. As we showed in Section r3.2l the spectrum has two main components, 
a continuum and a discrete component. Because of the Doppler-broadening the lines 
overlap, the result being a "continuum" with some little bumps in the position of the 
strongest lines. The line spectrum dominates only on timescales shorter than nuclear 
destructions one; over longer timescales the dominant part of the spectrum comes 
from the proton-neutron plasma radiation, which is mainly p-n Bremsstrahlung and 
p + n — » D + 7 capture. 

The destiny of a long-lived (i.e., t > r„), constant density and temperature 
plasma, is a purely protonic composition (cf. left panel of Fig. [7J . This is due to the 
/3-decay of free neutrons into protons, and as a result no other element will remain 
in equilibrium. Ref. [3| suggested a "neutron evaporation" mechanism through which 
a fraction of the neutrons escape the plasma, so that high-temperature, low-density 
plasmas may be astrophysical neutron sources. The neutrons might be captured 
by surrounding (relatively) cold(er) material such as hydrogen and form deuterium, 
emitting 2.22 MeV 7-rays. The observation of this spectral line would serve to reveal 
the existence of the nearby neutron source. 

Plasmas which for some reason are cooled to such an extent that reactions are 
frozen at the right moment in time, could produce a higher abundance of interme- 
diate elements (and therefore light elements) as shown in the right panel of Fig. 
[7j We conclude by suggesting that such a scenario could be physically plausible in 
situations where the plasma expands and cools, for example when being crossed 
by shock waves formed during a supernova explosion. Since supernova explosions 
are ubiquitous in the universe, this mechanism could, in principle, change cosmic 
abundances. 
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Appendix A. Evaluations 

Appendix A.l. Justification of the assumptions 

Here we intend to justify numerically the assumptions described in Section [2] 

First of all, we have restricted ourself to a domain of low densities and high 
temperatures, i.e. n ~ 10 18 cm -3 and T > 1 McV. As we have already mentioned 
in Section [TJ under such conditions break-up processes are dominant. Therefore, 
all heavy nuclei will be destroyed "instantly" , resulting in a proton- neutron plasma 
with some traces of light nuclei such as D, T, 3 He, 4 He. 

Secondly, the astrophysical environments where these type of plasmas are sup- 
posed to exist are such that we are allowed to start our calculations with an initial 
plasma composition similar to solar or cosmic composition of elements!^ 

The fact that all heavier nuclei are "instantaneously" destroyed, and that the 
initial composition is solar/cosmic (i.e. abundances of nuclei heavier than He are 
at least three orders of magnitude below) , suggest that heavier nuclei than 4 He will 
never become abundant. Below we attempt to quantify the importance of heavy 
nuclei interactions. 

First of all, we want to illustrate the claim that break-up processes are the most 
important processes for T > 1 MeV plasmas. In |Appendix B| for some important 
nuclei, are compared total nonelastic with inelastic cross-sections, and their corre- 
sponding reaction-rates. In that work are also compared the total nonelastic and 
capture reaction cross sections together with the corresponding rates. The plots 
show that the rates of break-up processes (the difference between total nonelastic 
and inelastic) is bigger than excitation (inelastic) and capture reactions for T > 1 



Secondly, by starting with a plasma composition similar to solar/cosmic chem- 
ical abundance, it is clear that the most abundant elements are hydrogen and he- 
lium, followed three to four orders of magnitude below by Oxygen, Neon, Nitrogen, 
Carbon, etc. Therefore, the major contribution will come from proton and alpha 
interactions due to their higher abundances. 

To quantify it, let us look at the destruction timescales e.g. for the 16 nucleus, 
as the third most abundant element. Let us consider reactions of proton, alpha and 
16 projectiles with 16 targets. The destruction reaction rates will be given by: 



MeV. 



^(p+ 16 0) = < c(p+ 16 o) v > n p 

F( a + le O) = < cr (a+ le O) V > fl a 

F(16 0+ 16 ) = < (7(16 0+ 16 ) V > 7J16Q. 



(A.l) 
(A.2) 
(A.3) 
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From Ref . l34l we get the ratios: n a /n p w 8 x 10 -2 , nie /n p ~ 5.4 x 10~ 4 . 

To estimate the order of magnitude value of the reaction rate for the above 
reactions we simply considered constant cross sections with value equal to the peak 
value. From measurements in Ref. [3"51 - [3T1 we get these peak values: <j( p+ i6 ) ~ 0.5 
b (total nonelastic), <J( a +w ) ~ 0.9 b (total nonelastic), and <J(i&o+ 16 o) ~ 1 b (it is 
of the order of the sum of all reaction cross sections) . 

The reaction rate for constant cross section is < a v > = a < v > where < v > 
is the mean relative velocity, < v >~ /i -1 / 2 , where /i is the reduced mass. Hence, 
the destruction reaction rates ratios are: 



r 



r 



(p+ le O) C(p+ ls O) I/ M(a, le O) 



^ ~ io- 1 



(16 0+ 16 ) (T(16 0+ 16 ) / ^( pj 16 ) 



r 



(p+ieo) ^tp+ieo) V M( 16 0, 18 0) V n v 



10^ 4 



The last estimate shows that neglecting interactions with projectiles different 
from protons and alphas will cause an error 

r(16 ° +16Q) < 0.1% 

r ( P + i6 o) + r( Q+ i6 ) + r(i6 0+ i6 0) 

Then, with an accuracy better than 0.1% we may neglect all nuclear reactions 
except those that involve protons and/or alphas. One must note here that due to 
the break-up processes, the number of neutrons increases very fast and they be- 
come very abundant. Therefore neutron interactions become important soon. After 
taking all these into account, we conclude that the most important nuclear inter- 
actions in high-temperature, low-density plasmas, are those reactions which involve 
protons, neutrons and alphas as projectiles. The interactions of heavy nuclei will 
never become important. 



Appendix A. 2. Photodisintegration 

As discussed in | Appendix A.i] the most important projectiles are protons, neutrons 
and alpha particles. Plasmas produce photons through ion-nuclear reactions and 
electron Bremsstrahlung (see Section. [5]). The photon number density depends on 
the optical depth, i.e. electron density and the plasma's size. The photon number 
density per E 1 energy interval (dn^/dE^), is given by: 

dn 7 dh~ R 

^^ x 7 xmax(1 ' T) ' (A - 4) 

where dhj/dEj is the photon production rate per E 1 energy interval, R is the size 
of the plasma, r is the optical depth and c is the speed of light. The term max(l, r) 
takes into account the photon's time-escape delay due to propagation through an 
optically thick plasma (r > 1). 
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The optical depth for a homogeneous plasma with number density n and size R 
is given simply by t{E 1 ) — n R a c (E 7 ), where a c (E 7 ) is the Compton scattering 
cross-section (Klein- Nishina). 

The photodisintegration reaction-rate will be given by: 



T( 7 ,x) = / dE y ^a h . x) (E 7 )c (A.5) 



where a^ ltX ^(E 7 ) is the photodisintegration cross section for a given nuclei. The 
cross sections for the important nuclei such as D, 12 C, 16 0, 56 Fe, etc, generally 
show a threshold E 1 > 3 MeV, and have a maximum around 7 — 20 MeV. 
By including all of the above, one obtains: 

r (7 , x) = R JdE^ ( dT ^f r) ) x a M (E y ) x max(l,T(£ 7 )) (A.6) 

Note that the photodisintegration reaction rate in formula IA.61 includes two 
cases, one without optical depth and one with optical depth. The case without 
optical depth scales as n 2 R (n e is the electron number density), whilst the case 
with optical depth scales as R 2 . On the other hand ion- nuclear reaction rates 
scale as n 2 (rii is the ion number density). Thus, for an optically thin plasma, 
photodisintegration effect is mostly influenced by the size R because n e ~ m. 

In high temperature (Tj > 1 MeV), low density, optically thin (r < 1) plasma, 
the electrons are heated due to Coulomb exchange with ions and are cooled due 
to Bremsstrahlung, which has a higher rate than heating. In the extreme case of 
electrons getting heated to high temperatures (T e > 1 MeV), pair-production pro- 
cesses become non-negligible, and thus, the plasma will cool down (see Ref. [52] and 
Ref. [53]). 

The electron thermal relativistic Bremsstrahlung (e-e and e-ion) emissivity is 
given in Ref. 1511 The p-n capture and Bremsstrahlung thermal emissivities are 
calculated as described in Section. 13.21 

We consider here an uniform, electrically neutral, spherically symmetric, opti- 
cally thin plasma, with a size R ~ 10 6 cm and number density n ~ 10 18 cm~ 3 . 

Two important nuclei are considered, the 16 O as the third most abundant cos- 
mic/solar element after Hydrogen and Helium, and Deuterium as the most abundant 
light element obtained after heavier nuclei are destroyed (see figures of Section. [3]) . 
By using the formula fA.61 we calculate the reaction rates for 16 0(7, X) and D(7, n)p 
photodisintegration reactions. For the sake of comparison, we calculated separately 
the photodisintegration reaction rates of photons emitted by electrons and ions, 
and we include the 16 0(p, X) proton break-up. The results are shown in Fig. 8. 
We notice that reaction rates coming from electron Bremsstrahlung photons are a 
function of electron temperature, whereas the reaction rates coming from nuclei are 
a function of ion temperature. 
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t 1 1 1 1 1 1 1 r 




T (MeV) 

Fig. 8. The photodisintegration reaction rate for ls O and D nuclei as a function of temperature, 
for a plasma with size R = 10 6 cm, and n = 10 18 cm -3 . The full line ( ls O) and the dotted curve 
(D) are the photodisintegration reaction rates of photons comming from electron Bremsstrahlung 
as a function of electron temperature (T e ). The dash-dotted ( 16 0) and the full-dotted (D) lines 
are the photodisintegration reaction rates of photons coming from nuclei as a function of ion 
temperature (T{). For comparison the dashed curve is the proton destruction reaction rate of 16 O 
as a function of ion temperature (T;). 

Fig. 8 reveals that for a spherically symmetric, optically thin plasma with size 
R = 10 6 cm, and number density n = 10 18 cm -3 , and for a solar/cosmic initial 
abundance, the photodisintegration rate is negligible. Even when the electron tem- 
perature becomes comparable with the ion temperature, the photodisintegration 
reaction rate is 3-4 orders of magnitude below other nuclear reaction rates, which is 
clearly seen by comparing 16 0(7,X) and 16 0(p, X). In the case when the electron 
temperature is lower than the ion temperature, and there is a significant amount 
of neutrons in the plasma (here we took, for simplicity, n n — n p — 10 18 cm" 3 ), 
the most significant amount of 7-rays are produced from p — n interactions. We see 
from Fig. 8 that the contribution of photodisintegration due to 7-rays coming from 
p — n interactions is negligible. We note that nuclear prompt 7-ray lines have no 
significant effect in photodisintegration. The same is valid for Deuterium, except if 
it is in nuclear statistical equilibrium (NSE) as it is in Section. [3] For this case, the 
equilibrium abundance will strongly depend on electron temperature. If the electron 
temperature is T e < 1 MeV then photodisintegration due to photons emitted by 
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electrons does not play any role, whereas the photodisintegration by photons com- 
ing from p — n interactions could affect the NSE equilibrium for timescales longer 
than 10 4 s. 

If the size of the plasma and/or electron number density is smaller, photodis- 
integration will be negligible. However, if the size of the plasma is R > 10 s cm, 
and/or its number density is n > 10 18 cm -3 , or/and the electron temperature is 
very high (T e > 1 MeV), then photodisntegration must be taken into account. 



Appendix A. 3. The relevant projectile energy interval 

Currently we do not have a complete nuclear theory which allows to calculate pre- 
cisely the a\j for any requested energy interval. Therefore, a good approximation 
must be used here. 

We are interested in plasma reactions in the temperature interval 0.1 — 20 MeV. 
For any given temperature, there exists an energy interval or window, such that, 
integration over it, brings the major contribution to the Maxwellian average (for 
that particular temperature). 



,00 , £l (T) 




p 00 


/ = / + 
10 JO 


7 ei (T) _ 


+ / 



C iT) + 1 



^2(T) _ 



p 



0.1% 



(T) 



Working within the temperature range of 0.1 — 20 MeV, we want to obtain the 
range of the integration [E\,E2] such that the accuracy of the Maxwellian average 
is greater than a desired value (e). The energy interval [E±,E2], will contain any 
energy window [ei(T), ^{T)] for any temperatures in 0.1 — 20 MeV. It is expected 
that the lowest temperature (Ti) will affect the value the lowest energy E\, and the 
highest temperature (T2) that of the energy E%. This is true due to the distribution 
functions of these two extreme temperatures. 

Bearing in mind the above argument, we can estimate E\ and E 2 from T\ and T 2 , 
separately. For simplicity, let us consider an extreme case of constant cross section 
with a value equal to the peak of the actual cross section. Since the energy window 
upper limit (£2) for T\ and the energy window lower limit (ei) for T2 lies inside 
[Ei, E2] we can make the following approximations to the integral in Eq.([T]): 



/„ 1 dEE exp(-E/Ti) 
J™ dEE exp(-E/Ti) 

J~ dEE exp(-E/T 2 ) 
J E2 dEE exp(-E/T 2 ) 
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Integrating the above equations we obtain the following algebraical equations 
that allow to find the root: 

(1 + ^ )cxp( ^ ) = TT7' 

If we take Tl = 0.1 MeV, T 2 = 20 MeV, and e = 0.1%, the solutions are 

E 1 w 0.015 MeV E 2 ~ 133 MeV 

For less accurate results, the interval [E\, E2] will be narrowed. For the capture 
reactions such as (n, 7) , where the cross sections at low energies show a dependence 
a ~ 1/v, it will be a more complicated function but the lower valid energy holds 
for e = 1%. 

The last estimation informs us numerically of the energy interval beyond which 
the contribution is negligible. This interval is E G [0.01, 150] MeV for accuracy 
< 1%. The better we know the cross sections in this interval, the more accurate the 
values of the reaction rates will be. 

Appendix B. Cross sections and reaction rates for some typical 
nuclei 

Appendix B.l. Total nonelastic and inelastic cross sections with 
their corresponding rates 
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Appendix B.2. Comparison between inelastic and capture cross 
sections for some specific elements 
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Appendix B.3. Production of j-ray lines for some elements 
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